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' We present a self-consistent field theory (SCFT) of semiflexible (wormlike) diblock copolymers, 

, each consisting of a rigid and a flexible part. The segments of the polymers are otherwise identical, 

■ in particular with regard to their interactions, which are taken to be of an Onsager excluded- 
' volume type. The theory is developed in a general three-dimensional form, as well as in a simpler 

one-dimensional version. Using the latter, we demonstrate that the theory predicts the formation 
j^'. of a partial-bilayer smectic-A phase in this system, as shown by profiles of the local density and 

orientational distribution functions. The phase diagram of the system, which includes the isotropic 
and nematic phases, is obtained in terms of the mean density and rigid-rod fraction of each molecule. 
The nematic-smectic transition is found to be second order. Since the smectic phase is induced solely 
' by the difference in the rigidities, the onset of smectic ordering is shown to be an entropic effect 

and therefore does not have to rely on additional Flory-Huggins-type repulsive interactions between 
unlike chain segments. These findings are compared with other recent SCFT studies of similar 

■ copolymer models and with computer simulations of several molecular models. 
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I. INTRODUCTION 



Liquid crystals, being of major importance both on fundamental grounds and for industrial applications, have long 
■ been of interest to researchers 0, |[. Another class of materials of considerable interest are melts of block copolymers, 
Ph ' which have been shown in recent decades to exhibit a multitude of phases of varying complexity . This paper will be 
O . concerned with liquid-crystalline copolymers, i.e., copolymers which exhibit phases commonly associated with liquid 
crystals. 

For the study of nonuniform polymer melts, a particular type of mean-field density-functional approach traditionally 
known as self-consistent field theory (SCFT) has proven to be a very powerful tool [§, |[ |[ @, || ■ Aside from obvious 
determinants like the composition of the melt being studied, the models considered by SCFT approaches differ mainly 
in two regards: the flexibility of the polymer chain and the interactions between chain segments. Most commonly, 
polymers have been modelled as perfectly flexible Gaussian chains, which assumes that there is no energy penalty for 
local bending. Alternatively, the "wormlike" chain model |l^ does introduce such a bending penalty and is thus the 
appropriate model for semiflexible polymers. Concerning the interactions, most previous investigators have invoked 
CSj Flory-Huggins-type repulsive interactions between unlike chain segments in order to induce microphase separation 
: and hence the formation of various mesoscopic phases |ll|, |l^ . 

The aim of the present paper is to develop a modified SCFT suited to examining the formation of liquid-crystalline 
phases in melts of diblock copolymers each consisting of a rigid ("rod") and a flexible ("coil") part. Both parts 
are modelled as wormlike chains, but are characterized by different rigidities. For the interactions between any two 
chain segments, in this work we do not distinguish between the two parts, i.e., we assume that the same type of 
^ ' interaction applies to any pair of chain segments. We take this interaction to be the limit of the Onsager excluded- 
Q , volume interaction for thin rods, which favors local alignment of the chains. In the case of homogeneous chains, 
O ■ characterized throughout by the same rigidity, the present model reduces to that applied by Chen et al. Q to 
^ ' the study of nematic ordering and the isotropic-nematic interface of semiflexible polymers. Here we show that the 
generalized model with different rigidities for two parts of a polymer is able to account for the formation of lamellar 
smectic-A phases. With this minimal model, the only possible sources for the formation of such phases are entropic 
^ , effects. Recently, both smectic monolayer and bilayer phases (traditionally denoted Ai and A2) have been found 
■ - ' in an SCFT-based approach which includes Flory-Huggins interactions and treats the rigid sections of the polymers 
as perfectly rigid and perfectly aligned in the same orientation [p^ [l^ . In contrast, here we find only one type of 
smectic-A phase, which would most accurately be classified as a partial bilayer phase (denoted Ad). In addition, as 
in Refs. |l^, [ij], our model accounts for the disordered isotropic phase. 

Liquids of fairly short rod-coil molecules have been examined recently by computer simulations |l8|, |l^, |2C , 21 , 22j . 



Generally, these computer studies have shown that the addition of flexible segments to otherwise rigid molecules 
stabilizes the smectic-A phase with respect to the nematic phase, consistent with experimental results for thermotropic 
non-polymeric liquid crystals 1^. This behavior is also shown by the present model, although, in contrast to some 
studies [ p8| , [lS| , |2^, we do not find that the nematic phase is entirely suppressed with respect to the isotropic and 
smectic phases. Our finding that the only stable smectic-A of the present model is the partial bilayer Ad phase is in 
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agreement with a density-functional treatment by Holy st l 23[ of rigid "nail-shaped" molecules, and, somewhat more 
tentatively, with Monte Carlo studies by Mazars et al. [p_9| of a rod-coil model. 

Experimental studies by J. T. Chen et al. demonstrated both monolayer and bilayer smectic ordering in 

rod-coil diblock copolymers consisting of a highly rigid polyhexyl-isocyanate block joined to a flexible polystyrene 
coil. However, these results were complicated by the simultaneous occurrence of tilted (or smectic-C) ordering and 
crystallization of the rigid blocks. These two effects are not considered here, but will be examined in future work. 

This paper is organized as follows. In Section II, a detailed description of the theory is given. In Section III, 
numerical results are presented, showing the local density and orientational order-parameter profiles in the smectic 
phase as well as the phase diagram of the model in terms of mean density and rod fraction. We conclude with a 
summary and discussion in Section IV. 



II. THEORY 



General 3-dimensional theory 



We consider a monodisperse melt of n rod-coil diblock copolymers of total contour length L, polymerization index 
N, and fixed segment length a such that L = Na, occupying a total volume V. A fraction a of the total contour 
length of each copolymer is occupied by relatively rigid (rod) segments, and the remaining fraction by more flexible 
(coil) segments. The total average monomer density po = Nn/V = Np is fixed, whereas we shall not require the 
melt to fulfill a local incompressibility condition. In accordance with the wormlike chain model for semiflexible chains 
1^, |l^, which is applied here to both sections of the chains, polymers will be treated as space curves 

ri{t) characterized by dimensionless unit tangent vectors Ui{t). Hence, microscopic contour- averaged density operators 
may be defined as 



-1 ^ pa 



(t)rigid{r, u) = 7 7 ^ / dti 6{r- ri{ti)) S{u- Ui{ti)) 

(^/iex(r,u) = - V / dU 6 {r - r^{U)) S (u - u,{U)) , (1) 



which satisfy the normalization conditions 



dr du 4)rigidir,u) = aV, 

drdu(^/ie.(r,u) = (l-a)y. (2) 

Since these partial densities are only determined a posteriori and otherwise do not enter the calculations, we shall 
instead use the total density, 

(f){r, U) = (f)rigid{r, U) + (f)flexir, u). (3) 

The interaction potential between any two segments, i.e., rod-rod, rod-coil, or coil-coil, is taken to be the Onsager 
excluded- volume interaction, which in the limit of very thin polymers (i,a >> diameter D) reduces to [14 

D(ri,Ui;r2,U2) = L'^D S{ri - r2) |ui x U2I. (4) 

The partition function in the canonical ensemble has the form 

Z = JVn{-}exp(^-pC J dr J duj du' ${r,u) ^{r,u') \uxu'\^ , (5) 

where 



C = L^Dp 



(6) 
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is proportional to the average polymer number density p. In Eq.(|^), 

j !?„{.}= 1 J] j I?{r„uj7'{r„u,[0,l]}, 



(7) 



where 



P{ri,Ui[si,S2]} 



S2 r f' 1 1 f 

\{5[n,{tf-l\5 r,{t) - T,{si) ~ L ds u,(s) exp -— / 



dt' K(t') 



du,: 



(8) 



is the statistical weight of a given path and Hi{t) is a dimensionless bending modulus. According to the model 
considered here, 



n{t) = Krigid , < t < a 
= Kflex , a <t <1 . 



(9) 



As is a standard procedure in self-consistent field theories, we now multiply the partition function by 1 = 
/ 1? {0(r, u)} (5((/)(r, u) — (f>{r,u)), which allows us to replace the operator <j) with the function (f>. Then, us- 
ing the exponential representation of the delta function, the partition function can be rewritten as Z oc 
/ V{W} f V{(j)}exp{-J'[W,(f\) with 



C / dr / du / du' (/)(r,u)0(r,u')|u x u'| 



dr / du W{r,u)(f>{r,u) 



In 



where 



Q = J Vi{-} exp dtW{r{t),u 



(10) 



(11) 



is the single-chain partition function. The function W{r, u) is identified with the potential energy, or field, generated 
by all polymers in the system as seen by a "test" polymer. 

Taking the saddle point of the free energy (|l^ with respect to (/)(r, u) and M^(r, u), we obtain the mean-field 
equations: 



W{r,u) = 2Cy" du' 0(r,u')|u X u'l, 
V SQ 



(r,u) 



Q SW{r,u)' 



(12a) 
(12b) 



In this approximation, the function 0(r, u) equals the statistical average < 0(r, u) > of the microscopic density. In 
order to solve these equations for the density, we need to express (5[VF] in terms of the end-segment distribution 
function defined by 



q{r,u,t)= Jv{r,,u,}r{v„u,;[0,t]}S{v~r,{t))S{u-u,{t))exp - ds W (r,(s), u,(s)) 



(13) 



The function q^(r, u,t) is defined as the end-segment distribution function starting from the opposite end of the 
polymer. Hence, 



Q ^ J du J dr q{r, u, i)gt(r,u,t). 



(14) 



where the contour variable t is arbitrary. With the above definitions of Q, q, and g^, Eq.( 12b ) yields for the local 
density: 



V . 
'")^7^/ dt (7(r,u,t)(7T(r,u,i). 
^ Jo 



(15) 
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The average rigid and flexible densities, (/)rigi(i(r, u) and (f>fiej;{r,u), are given by expressions analogous to (|5|) on 
replacing the limits of integration over t as in Eq.(|l|). The end-segment distribution functions, or propagators, satisfy 
diffusion-like equations: 



d_ 
di 

d_ 

di 



q(r,u,<) = 
g'f(r,u,t) = 



-iu • Vr + 

-LU • Vr - 



1 



1 



g(r,u,i), 
gt(r,u, i). 



(16) 



with initial conditions (7(r, u, 0) = 1 and (/^(r, u, 1) = 1. Here we have defined the rigidity parameter ^(i) = K{t)/N 
depending on the chain contour variable t: ^{t) equals the persistence length of the corresponding chain section in 
units of the total contour length L This is where the difference between the rigid and flexible parts of the 

copolymer is accounted for. 

To proceed with further analysis of the mean-field equations, we represent the orientational (u) dependencies of 
the functions cj), W, q and using spherical-harmonic series: 

'/'(r.u) = ^0/„(r)yi^„(u), 



q{r,u,t) = ^g/™(r,t)y,,™(u). 



l,m 



gt(r,u,t) = J2qljr,t)Yi,m{n). 

l,m 

Since these are all real functions, the expansion coefficients must obey the following conditions: 

etc. Next, we expand the kernel |u x u'| by use of the addition theorem for spherical harmonics fl^ , 



\nxu'\^J2 



47r 
21 + 1 



di y,,™(u)r,* (u') 



with 



di — Q , I odd, 

7r(4fc-hl)(2/c)!(2fc-2)! 



^0 = 1 



d2k 



k= 1,2,3,.. 



Inserting these formulae into the free energy (^0|), the latter can be expressed as 

r AttC 



^^pY. dr[^^d,|0,,,„(r)|2 - Re (W^,,„,(r)</>r,„,(r)) 



- In 



where "Re" denotes real part. The mean-field equations (12a) and ( p^ ) now read 

Stt 



21 + 1 



diC (jji^mir) , 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



</'i,m(r) 



V 



V .nv L" .m" 



^/ E 9/^,„'(r,i)gr,™"(r,i)i 

^ I'm' I" ,ni" 



' i2l" + i)i2l' + i) 

An{2l + 1) ^™",™',m > 



(23) 



5 



with 



(24) 



The C^;, ^, ^ are Clebsch-Gordan coefRcients, and we have used a result for the integral of three spherical harmonics 
|p5| . The corresponding projections 4>rigid,i,m and 4>fiex,i,m are obtained likewise by changing the limits of the line 
integral in Eq.(23). In terms of the projections qi^mi^^t), the diffusion-like equation ( p^ yields the following coupled 
set of equations: 



2/' + 1 



^ E V 2/ + 1 

n' 

(f^l,l',l _ 



^1,1', if }_f^l.l'.l \ d 

•-^0,0,0 \^ y-'l,m',m,~^ ^-l,m',mj 



1 



c 



1,1' I 



V2 

d_ 

dy 



dx 



E 



lil + l)qi,m{r,t) 



I' I'- 



' (21" + l){2l' + I) „irt)&"^^c''^^ 



with initial conditions 



go,o(r,0) = V47r , 
gi,m(r, 0) = 0, otherwise. 



A similar equation applies to q^ (r, u, t) . 



(25) 



(26) 



1-dimensional theory 

For simple applications of the general theory presented above, we now specialize to situations where the densities 
vary in only one spatial dimension, which for convenience is chosen to be the z direction. Furthermore, we will 
restrict analysis to phases that exhibit no azimuthal orientation dependence about the z-axis, thereby excluding the 
possibility of smectic-C phases |ll6| , Then the only nonzero projections of any angular functions are those with 
TO = 0, so that we will subsequently drop the to indices, denoting (pi = <j)i.rn=n, etc. Another consequence is that now 
all spherical- harmonic expansion coefficients are real. The free energy Eq.(|2l|) becomes 



jr = pAY^ / dz 



I 



21 + 1 



di(t>f{z)-Wiiz)Uz) 



In- 



,1 



(27) 



where A is the cross-sectional area of the system in the x and y directions. The mean-field equations (|2^) and (|2. 
are: 



Wi{z) 



21 + 1 



diC 4>i{z), 



t)qi"{z,t)A 



\2l" + 1){21' + 1) 



47r(2/ + 1) 



with 



Q = / q\{z,t)qi{z,t) 



(28a) 
(28b) 

(29) 



Finally, the diffusion-like equation ( p5| ) becomes 
d 



dt 



- E 



2m 



l{l + l)qiiz,t) 



' {21" + l)i2l' + 1) f,"j'j 



I' I" 



47r(2Z + 1) 



^(cSo') W,{z)qi„{z,t) 



(30) 
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with initial conditions 



90 (^,0) 



/47r , 
J > 0. 



(31) 



In the following, we shall truncate the interaction expansion (|l^) after I = 2, (equivalent to a "Maier-Saupe" 
interaction), retaining only the terms with coefficients do = f a-nd d2 = — ■§§ • 

We shall now briefly discuss the computationa l me th ods used in solving the theory. The fields and densities 
are determined self-consistently according to Eqs.(28a), ( ^8h| ) and (|o|) using a fixed-point iteration algorithm with 
variable mixing parameters for successive iterations. An iteration consists of: (a) given the functions Wi{z), solving 
Eq. ( |30| ) and its counterpart for qj ; (b) calculating the set of functions (jiiiz) from Eq. (28b); (c) obtaining a new set 
of functions Wl{z) from Eq. (28a); (d) mixing Wl{z) and Wi{z) according to a fixed-point algorithm, which yields a 
new Wi{z) for the next iteration. 

Concerning step (a), solutions of the diffusion-like equation ( |30| ) were discretized in time and space according to 
a Forward Time Centered Space (FTCS) scheme [ p6[ , which is explicit and straightforward to implement. This is 
a reliable discretization scheme for solving partial differential equations although its error is only first order in the 
contour discretization dt. For the corresponding diffusion equation in the case of Gaussian chains, a number of more 
sophisticated schemes such as Crank-Nicholson and DuFort-Frankel |2^ are available to improve on stability and 
accuracy (the error becomes second order in dt in both cases), and thus to lower the contour resolution needed for a 
given desired level of accuracy, a key factor in computing times. Due to the coupled nature of Eq.(^0|), however, it 
seems daunting to apply similarly efficient discretization schemes to the wormlike chain model. 



Co ncerning step (d), the fixed-point iteration algorithm used in finding self-consistent solutions of Eqs.(28a) and 
( ^8b| ) is a tried and true method. Here viable alternatives do exist, however, first and foremost being Newton- 
Raphson-type algorithms like Broyden's method (e.g. [^). In the present work, we have restricted ourselves to using 
the fixed-point algorithm, since our primary aim was to develop the framework presented above and illustrate it with 
only a few, demonstrative results. 

In the present work, calculations were performed on a 1-dimensional grid with periodic boundary conditions, a 
spatial discretization of dz = 0.02, and a contour discretization of dt = 1/1500. The spherical-harmonic expansions of 
the propagators and densities were truncated after I = 12. With these parameters, J-'/V can be determined to within 
a numerical inaccuracy of less than 1%, requiring up to 500 iterations, or 3-4 hours on a Pentium II processor. Note 
that although the interaction kernel (|l^) is truncated after / = 2, we cannot truncate the propagators at the same 
value of I, si nce t he cou pling i n (3C| ) renders these higher-order projections nonzero. 

Once Eqs.(28a), ( ^8b| ) and (^0) have been solved self-consistently, the free energy can be rewritten as 



J' = -inCpAY^ 



21 + 1 



dz (j)f{z) — In ■ 



(32) 



For determining coexistence regions of the first-order isotropic-nematic transition, it is necessary to perform double- 
tangent constructions on curves of the free energy per volume J-'/V. Omitting linear terms in C (which yield constants 
when the derivative is taken), we obtain 



T TC 



-4.C^ ( ^ 



^ 21 

I 



- y dz ^f{z)--C\n^ + ClnC, 



(33) 



where we have used Stirling's approximation for the factorial. In a smectic phase, all projections of the local densities 
and propagators are taken to be periodic functions of z with period (or layer spacing) D. All integrations over z can 
then be evaluated as 



A 
V 



dz 



dz 



(34) 



The location of the second-order nematic-smectic transition (cf. section III) is determined by the behavior of the 
smectic order parameter defined as: 



J dz ([V^(j)rigid,Q - a]^ + [V^tpflexfi - (1 - a)]^) 



1/2 



(35) 



The parameter Ogm vanishes in the isotropic and nematic phases, while it is nonzero in the smectic phase. For a 
"perfect" smectic phase in which the profiles \/^(t>rigid,o and V^4>fiex.o have rectangular shapes of widths aD and 
(1 — oi)D, respectively, Ogm has the value ^J2a{\ — a). The equilibrium period is that value of D which minimizes 
the free energy per volume. 



7 



III. RESULTS 



The results presented here apply to a fluid of polymers with S^rigid ~ 10, and ^/jei = 0.1. These values more or less 
represent the feasible bounds on ^ for performing numerical calculations. To begin with, we shall examine a system 
with a = 2/3. 

Fig. d shows the rigid and flexible density profiles (i.e., zeroth-order projections 0^.;=o(-2), P G {rigid, flex}) as well 
as the total density profile 4>o{z) in a smectic configuration ai C — 20, just above the nematic-smectic transition. We 
see a clearly defined region of predominantly rigid chain segments and a less pronounced region where the flexible 
parts are dominant. The individual density variations are nearly pure sine functions, as expected close to a second- 
order transition. The total density in this lamellar structure exhibits maxima in the rigid regions and minima in the 
flexible regions. This indicates that the rigid segments pack more efficiently than fiexible ones due to their greater 
susceptibility to local orientational alignment. 

More insight can be gained from the orientational projections of the densities, (j)i3^i{z). Here we express these in 
terms of the z-dependent order parameters defined as, for each degree /, 



"^''^'^^ ^ WTl M^) ' ^''^ 

equivalent to the average of the V"^ Legendre polynomial P;(cos0) per segment oi species p, where 9 is the angle between 
a segment axis and the z axis. The lowest-order functions Pi.p^z) and P2,ii{z) corresponding to the configuration of 
Fig. 1^ are shown in Figs._^ and |. The second-order distribution P2,/3 characterizes the overall degree of orientational 
order of species /3, while Pi^/j indicates the average spatial direction in which segments fi are oriented. The behavior of 
the functions Pi^^(z) shows that both the rigid and flexible regions are divided in the middle into domains of positive 
and negative orientation along z, in the manner of a bilayer. We see that both P2, rigid and P2jiex exhibit maxima 
(minima) in the regions of high (low) rod density (cf. Fig. |^). The minimum in the flexible coil distribution P2jiex is 
more pronounced, while its maximal region is a slightly modulated plateau within the domain of high rigid density, 
suggesting that coil segments are oriented by interactions with neighboring rigid segments. Overall, the coils exhibit 
a much weaker degree of orientational ordering than the rods. 

The configuration shown in Figs. ^ ^ and ^ has an optimal period (i.e., yielding the minimum free energy per 
volume) of D — 1.30L for the chosen parameters a = 2/3 and C — 2Q. Keeping the same a but increasing C to 30 
results in a decrease of the optimal period to 13 = 1.24. This is accompanied by stronger segregation of the rod and 
coil regions and less purely sinusoidal behavior of the densities, as shown in Fig. ^. The corresponding orientational 
order parameters are given in Figs. |5| and ^, and are qualitatively similar to those for C = 20. We have found that 
the optimal periods range from I.IL for a — 0.3 to 1.3L for a — 0.75 for values of C close to the N ~ Ad transition. 
Smectic phases with values of a outside this interval could not be examined due to numerical difficulties. Attempts 
to generate solutions corresponding to monolayer smectic structures having D < L always converged to a uniform 
(isotropic or nematic) phase. 

Fig. shows the phase diagram in the C — a plane for rod-coil polymers with frigid — 10 and £,fiex = 0.1. At 
low mean densities C, the steric interactions do not suffice to generate an ordered phase: the system is isotropic. 
Upon increasing C, we encounter a first-order isotropic-nematic (/ — N) transition. The two branches delineating 
the coexistence region in the diagram were determined by constructing tangents to the J- /V graphs for the isotropic 
and nematic solutions with a spline interpolation, which amounts to calculating the values of C at equal chemical 
potentials and pressures. The lower a, the higher is the value of C at which the I — N transition occurs: this 
transition is driven primarily by the tendency of the rigid-rod segments to align and thus requires higher densities if 
the proportion of the rods is lowered. For a = 1 we can compare our data with that of Chen who uses the full 
kernel ( p^ ) and obtains an isotropic density of Ciso = 4.18 and a nematic density of Cnem — 5.33. If we calculate 
this transition retaining the kernel up to Z = 12, we obtain good agreement: Ciso = 4.16 and Cnem — 5.29. With the 
truncated kernel {I < 2) used in the majority of this work, the coexistence region becomes more narrow: Ciso = 4.89 
and Cnem = 5.35. 

A second transition, from the nematic to the smectic (Ad) phase, occurs at higher values of C. This transition is 
second-order within numerical uncertainties, as is indicated by the behavior of the smectic order parameter Osm, Eq. 
(^). Fig. ^ shows that Osm starts to grow from zero at a well-defined critical value of the density C. There is a 
minimum of the N — Ad transition density at roughly a = 0.55. On moving toward both lower and higher values of 
a, the critical values of C increase sharply until numerical difficulties prevent us from extending the graph further. 
Obviously, a well-balanced proportion of rod-like vs. coil-like segments facilitates forming a smectic phase: the former 
stabilize the rod-dominated portion of the density profile (the nematic microdomain, as it were), and the latter, the 
coil-dominated portion. Deviations from the optimal a to higher values of a decrease the entropic advantage of the 
coils, which has to be compensated for by an increased density. Lower values of a, on the other hand, destabilize the 



8 



nematic microdomain ordering of the rods. This is why, for even lower values of a, the I — N and N — Ad phase 
boundaries approach each other until intersecting at a = 0.32. Below this value of a, the I— N transition is preempted 
by a first-order I — Ad transition. For a = 0.3, we find a coexistence region between Ciso = 29.3 and Csmec = 31.4. 
Unfortunately, numerical problems prevented us from examining even lower values of a: it thus remains an open 
question by how much the / — Ad transition deviates from the preempted I — N transition. 



IV. SUMMARY AND CONCLUSIONS 



We have presented a self-consistent field theory for semiflexible copolymers in a general three-dimensional, as well 
as a one-dimensional form without azimuthal orientational dependence. Using this framework, numerical calculations 
have been performed to investigate the occurrence of a smectic-A phase in systems of rod-coil copolymers where the 
two parts of each molecule differ only in their rigidities. 

A phase diagram was established for S,rigid = 10 and S^fiex = 0.1, which correspond to very stiff rods and very 
flexible coils. At low densities C, an isotropic phase is present, which becomes nematic for larger values of C via a 
first-order transition. At even larger values of C, a partial bilayer smectic-A phase forms, whose period increases with 
increasing rod fraction a. The latter effect, while small in the examined parameter range, indicates that the coils can 
overlap with other coil segments as well as with the rods, resulting in interdigitation of the layers. This feature is also 
suggested by the rather diffuse, sinusoidal variation of the local density profiles 0/3, o- The shape of the smectic phase 
boundary in Fig. exhibiting a minimum mean density C near a rod fraction a = 0.55, is plausible and consistent 
with experimental deductions [Q. (See also |^.) Preliminary calculations indicate that lowering the rod rigidity 
frigid shifts all transitions to higher densities C. 

In the present work, only anisotropic excluded-volume interactions between chain segments, assumed to be the 
same for the rigid and flexible portions of each molecule, are taken into account. To generate smectic phases, we 
have shown that it is sufficient to distinguish rigid and flexible parts of a polymer by means of their rigidities or 
reduced persistence lengths ^, and thus have demonstrated that lamellar ordering is a purely entropic phenomenon. 
This is consistent with the findings of computer simulations [0 |l^, |l], |2^ and a density-functional treatment 
p3[ of diblock models employing only anisotropic repulsive intermolecular forces. In contrast, the smectic ordering 
found in the SCFT treatments of Refs. ||l5| , [l6| is driven primarily by isotropic Flory-Huggins interactions between 
unlike chain segments, which may well be present in more realistic models. Unlike here, the latter works predicted 
the formation of both monolayer and bilayer smectic phases, as well as strong segregation between rigid and flexible 
domains, which we attribute both to the quite different nature of the interactions adopted and to the allowance of 
local compressibility in the present theory: as conjectured in Ref. [p^ , compressibility effects "could dramatically 
stabilize the bilayer phase" . 

We expect that in a more accurate treatment of excluded-volume effects, the density C at the N — Ad transition will 
not increase so sharply with increasing a as indicated in Fig. |^. Note that the Onsager excluded-volume interaction 
Eq. is taken to be spatially local, and hence by itself does not generate smectic phases in the absence of differing 
chain rigidities. A more realistic non-local treatment of the repulsive interactions between chain segments should 
generate a smectic-A phase even in the case of homogenous semiflexible polymers ^ ^ ^ . 

Future extensions of the theory presented here will attempt to account for the non-local interactions mentioned in 
the preceding paragraph as well as for smectic-C phases, crystallization of the rods, and non-lamellar morphologies, 
all of which are expected to occur at larger values of C and more extreme values of the rod fraction a. 
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FIG. 1: Density profiles of a smectic configuration with a = 2/3, frigid = 10, ^fux = 0.1, C = 20, and period 1.3L. The solid, 
dotted, and dashed lines correspond to -s/iTT^o, s/^4>rigid,o-i and \fiK(i>fiexfi-i respectively. 
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FIG. 2: Orientational order parameters Pi, rigid (dotted line) and P2,rigid (solid line) for the same case as in Fig. 1. 
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FIG. 3: Order parameters Pijux (dotted line) and P2jiex (solid line) for the same case as in Fig. 1. 





5: Orientational order parameters of rigid segments for the same case as in Fig. 
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FIG. 7: Phase diagram for rod-coil copolymers with frigid = 10, ^fUx = 0.1. The isotropic and nematic phases are separated 
by a narrow coexistence region bounded by the solid and dash-dotted curves. The boundary between smectic and nematic 
phases is second order, indicated by the dashed curve. 




FIG. 8: Smectic order parameter vs. C for a = 2/3. 



